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ABSTRACT 

The inviscid, linear, non-axisymmetric, temporal stability of the boundary layer associ- 
ated with the supersonic flow past axisymmetric bodies (with particular emphasis on long 
thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. 
The eigenvalue problem is computed in some detail for a particular Mach number of 3.8, 
revealing that the effect of curvature and the choice of wall conditions both have a significant 
effect on the stability of the flow. 

Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far down- 
stream solution are obtained for the stability analysis and compared with numerical results. 
Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall 
conditions, are presented. We find, in general, important differences exist between the wall 
temperature conditions imposed in this paper and the adiabatic wall conditions considered 
previously. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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§1. Introduction and Motivation 


The recent resurgence of interest in high speed flight vehicles has 
rekindled interest in supersonic and hypersonic flows. One important 
area of aerodynamic study is that of boundary- layer stability and the 
possible transition to turbulence. Since turbulent flows result in 
considerably greater skin frictions and heat transfer coefficients than do 
laminar flows, any method by which a boundary layer may be stabilized 
is worthy of investigation. 

Generally it is found in the case of supersonic boundary -layer flows 
that inviscid disturbances are more important (i.e. more unstable) than 
viscous disturbances. (Here, we characterise inviscid disturbances as being 
those with wavelengths comparable to the boundary-layer thickness, whilst 
viscous disturbances possess much longer wavelengths). This is in contrast 
to the situation encountered for many incompressible bounda ry- 1 aye r flows 
where viscous instabilities are generally dominant. The first authors 
to give any form of rigorous mathematical account of the stability of 
compressible boundary layers were Lees and Lin (1). Making use of 
a rational asymptotic approximation they determined that the quantity 

[ p* j (where u* denotes velocity tangential to the surface, 
dy* 3y 

y* the normal to the surface, and p* the fluid density) plays a key 
role, similar to that the quantity 9^u*/9y*2 has in incompressible 
theory. In particular at the point where the above expression is zero 
(y*=yj*) i termed the generalised inflexion point, then there may exist a 
neutral mode with wavespeed u*(yj*); neutral modes are classed as being 
''subsonic'', ’'sonic'' or ’’supersonic'' depending on how the freestream 
Mach number is related to the wavespeed (Mack (2)). If the neutral 
disturbance is subsonic then the mode decays in the far-field; supersonic, 
neutral disturbance modes exhibit an oscillatory behaviour in the far-field; 
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a sonic mode occurs at the crossover point between subsonic and supersonic 
cases. Mathematically, these classifications are directly related to the 
non-dimensional wavespeed c (defined in (2.10) below), and the 
free-stream Mach number M,,. For subsonic disturbances we have 

1 • A." c < '♦ SL • 

for sonic disturbances we have 

c ' '■ R, or c ’ PL • 

and for supersonic disturbances we have 

c •= '• K °' C>, *L 

Arguments relating to generalised inflexion points have no implications 
for supersonic neutral modes. 

Lees (3) considered the effect that wall cooling has on the 
stability of compressible boundary layers on the basis of asymptotic 
theory. He predicated that with sufficient wall cooling the boundary 
layer could be completely stabilized and presented a criterion whereby 
the ratio of wall temperature to the recovery temperature at which the 
critical Reynolds number becomes infinite, can be computed. Even though 
Lee's original work contained numerical errors, subsequent authors 
including Van Driest (4) and Dunn and Lin (5) showed that Lee's predictions 
appeared to be correct. 

Van Driest (4) calculated the cooling required to completely 
stabilize the boundary layer on a flat plate at supersonic speeds with 
zero pressure gradient. Whereas Lee's investigations were limited to low 
supersonic flows. Van Driest predicated that complete stabilization was 
achievable by wall cooling over a wide range of Mach numbers up to 
hypersonic flow. He found, however, that for Mach numbers greater than 9, 
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it is impossible to stabilize the boundary layer with any amount of 
cooling when a Prandtl number of 0.75 and the Sutherland 
vi scosi ty- temperature law are assumed. 

The above predictions of Lees (3) and Van Driest (4) were based 
on the asymptotic theory of two-dimensional disturbances. Dunn and Lin 
(5) extended this work to include three-dimensional disturbances and 
they indicated that cooling was indeed an effective method by which the 
boundary layer could be stabilized for moderate supersonic Mach numbers. 
Based on their asymptotic analysis, Dunn and Lin, however (wrongly) 
concluded that at supersonic free-stream Mach numbers the boundary layer 
can never be completely stabilized by cooling with respect to all three- 
dimensional disturbances. 

However the above asymptotic work is found to have its limitations. 

The theory above indicates the possible existence of only subsonic 
generalised modes for cooled wall conditions. In a series of papers Mack 
(2, 6, 7, 8, 9, 10, 11), using numerical techniques, demonstrated that 
for compressible boundary layers there in actual fact exist a large 
(probably infinite) number of unstable modes, including many with supersonic 
neutral points, which may not be predicated by the above asymptotic work. 

For the case of three-dimensional disturbances in Blasius-type boundary 
layers, Mack showed that the so-called ’’first mode’’ of instability 
undergoes stabilization with sufficient wall cooling, thus verifying the 
above mentioned asymptotic analysis. However, Mack found in the case of the 
’’ second mode’’ of instability (which turns out to be more important because 
of the larger growth rates involved) cooling the wall actually destabilizes 
this mode. According to recently presented numerical observations by Mack 
(2, 10) oblique disturbances can be completely stabilized by wall cooling 
for Mach numbers up to 5.8, although they do require a larger amount of 
cooling than the corresponding two-dimensional disturbance terms. 
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Recently, Duck (12), considering the effect of curvature on the 
inviscid axisymmctric stability of a compressible boundary layer, 
determined that it has a profoundly stabilizing effect. Considering 
axisymmctric, inviscid disturbances in the supersonic boundary layer formed 
on a thin straight circular cylinder with adiabatic wall conditions, he 
determined that the effect of curvature is to completely stabilize the 
'’first mode’’ of inviscid instability at a comparatively short distance 
down the axis of the cylinder and causes the ’’second mode' 1 of instability 
to undergo substantial maximum growth rate reductions at increasingly 
further downstream locations along the cylinder axis. All the numerical 
evidence though, is that for large distances downstream the ’’second mode’’ 
is still present, having a very small, but measurable growth rate. Duck 
also observered that the inclusion of curvature terms causes the 
generalised inflexion condition described above to be altered and he 
derived a modified (or ’’doubly generalised’’) inflexion condition, 
involving the radius of curvature. 

The work of Duck (12) was extended to include non-ax i symme t r i c 
disturbance terms by Duck and Shaw (13) (here-after refered to as DS) 
and they applied their techniques to a different axisymmetric configuration, 
namely that of a sharp cone. DS generalised the doubly generalised 
inflexion condition further, to give a condition for the existence of so 
called inviscid, neutral non-axi symme t r i c modes and termed it the 
’’triply generalised inflexion condition 1 '. 

In this paper we consider the combined effects that curvature 
and wall cooling (or heating) have on the stability of the compressible 
boundary layer associated with an axisymmetric body (principly a thin, 
straight ci rcular x:yl inder) . We are interested in both axisymmetric and 
non-axi symmetric, inviscid disturbances. The theory for the form of the 
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boundary layer and the resultant disturbance equations has already been 
presented by Duck (12) and in DS, to which we refer the reader to for 
details, but we shall begin this paper with a brief account of the relevant 
equations to the present problem. It should be noted that although the work 
in DS was for a sharp cone, the necessary equations for a circular cylinder 
for this problem are easily obtained by setting the slope parameter, A. 
in DS to zero. 


§2 Formulation of Problem 

The general layout of the problem is shown in Figure 1. The z* 
axis lies along the cylinder axis, r* denotes the radial coordinate, 
and 0 the azimuthal coordinate. We have a supersonic flow incident 
normal on the end- face of the cylinder (which has radius a ), where 
Moo represents the f ree- st ream Mach number , the free-stream veloci ty 

(in the axial direction) and p„, and T„ represents the free-stream 
density, first coefficient of viscosity and temperature respectively. We 
define the Reynolds number to be 

Re = U ~ f -fi a , (2.1) 

and this will be assumed to be large throughout. Subscript » denotes 
free-stream conditions. 

For the resultant boundary layer formed on the cylinder we make the 
following assumpt ions : - 

Firstly following previous authors we shall assume an absence 
of any shock wave in the basic flow. We feel justified in making this 
assumption, as the chosen ’thinness’ of the boundary layer is such that the 
shock wave will be sufficiently far away from the boundary layer to have 
negligible effect. Secondly we shall assume the linear Chapman viscosity 
law, where 

H* = C T\ (2.2) 

and C is assumed to be constant. Since the flow is one of high 
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Reynolds number, we make the steady boundary- 1 ayer approximation. A 
fundamental component of this work is the inclusion of curvature terms 
into the governing equations. To achieve this we assume that the body 
radius is generally of the same order as the boundary - 1 aye r thickness (a 
similar approach was adopted by Seban and Bond (14), Duck and Bodonyi 
(15), Duck (12) and in DS), except at the tip of the cylinder where we 
shall assume the boundary layer has zero thickness (as assumed by Seban and 
Bond (14)). The final assumption we shall make is that of a perfect gas. 

Making use of the above mentioned assumptions, in DS, we obtained the 
following system of non-dimensional, leading order equations for the basic 
flow 


Irlir 1 ] + T+ + I^[t 2 ] =0 ’ 

dP 


37 


= 0, 


V| 3? + VJ 1 r T fr ! ]• 

v .£ + v3£ = tW(|p ) 2 


(2.3) 

(2.4) 

(2.5) 


where Uoo(Re* J Cv j ,V 2 .vj) represent the velocity components in the 

(r*=a*r, 0, z* = Re a*C'^z) directions respectively, p*R*T*P 
♦ ^ 

(where R is the gas constant) and T M T represent the pressure and 

temperature terms respectively, y denotes the ratio of specific heats 

p C D * 

and a = — j— F is the Prandtl number, with K being the coefficient of 
K 

heat conductivity; in (2.3)- ( 2 . 6 ) the basic flow is assumed independent of 
0, and has no azimuthal velocity component (i.e. V 2 = 0 ). 

The boundary conditions for the problem are 
V 1 = v 3 = 0 on r=l , 

V 3 ■* 1 , T -* 1 as r-*», (2.7) 

together with a wall temperature condition which for heated or cooled walls 
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has the form 

T = T w on r=l . (28) 

The system of equations (2.3)-(2.6) was solved using a straight- 
forward Crank-Nicol son scheme, as described in DS, subject to the boundary 
conditions (2.7) and (2.8). 

To study the stability of this basic flow, we now investigate the 
effect of small amplitude disturbances. We shall assume that the 
disturbance wavelength is generally comparable to the boundary- 1 aye r 
thickness and therefore also of the body radius (0(a )), in which case 
the parallel flow approximation is completely vindicated. 

At a fixed z-station we express the flow parameters of velocity, 
pressure, temperature and density as the sum of a mean flow term plus a 
small, first order disturbance term, for example. 

vi = 5a(C <p( r) E + 0(8 2 ), 

V3 = U 00 [wo(r) + 8 V3(r)E] + 0(S 2 ), 

T* = t! [T 0 ( r) + 8 T(r)E] + 0(8 2 ). 

P* = pJI*T* [ l+8p( r )E] + 0(8 2 ) , (2-9) 

where 

E = exp[ ^ (z-ct) + in 0], (2.10) 

T 0 (r) = T (r.z) 

wq( r) = v 3 (r,z) (2.11) 

and 

c = zi, T1 = (r-1)/;. z = z*/a* , (2.12) 

and 8 is the perturbation parameter and as such is assumed to be 
dimini shingly small . 

In this paper we focus our attention on temporal stability for which 
the growth rate is acj, where a represents the spatial wavenumber 
and Cj is the imaginary part of the wavespeed. If c s > 0 we then 
have growing disturbances, if c j = 0 the disturbance is neutral and if 
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Cj < 0 the disturbance decays. 


Substitution of the flow parameters into the full system of equations 
of continuity, momenta, energy and state, discarding 0(5 2 ) terms and 

all but the largest terms in Re, we obtain a sixth order system, which 
after some algebra (as described in DS) can be reduced to the following 
system 


jg 2 ( w o- c ) 9 = . h— , 

T 0 yM 2 ® 


(2.13) 


~~ <J> = (w 0 -c)q>q - [Wq^ - , ( W Q - c } C ) <p 
ym 2 „ 1 i+nc 


whe re 


♦ = T 0 [1 + ~ - I - M^(w 0 -c) 2 . 

ct 2 (l+riC) 2 

The boundary conditions are 
9 = Pq = 0 on q = 0, 


and 


(2.14) 


(2.15) 


(2.16) 


i a 

<P - 9® K n (ii) 
p - p®M»i) 


as t| 


(2.17) 


where 

q = ± a[ 1 -M^( 1 - c ) 2 ] i (£ + q). (2.18) 

It should be noted that the argument of the modified Bessel function 

(i.e. appropriate sign in (2.18)) K n (q) is chosen to ensure boundedness 
in the farfield. 

The system of equations (2.13), (2.14) was solved using a Runge— 
Kutta scheme for the eigenvalue c. as described in DS, subject to the 
boundary condi t ions (2.16), (2.17). 
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§3 Numerical Results 

Most of this section will be devoted to the effect of wall-cooling 
on the stability of compressible boundary layers and its interaction with 
curvature, as these conditions are likely to be of interest in important 
applications. However, at the end of the section we shall give some heated 
wall results which exhibit some additional interesting physical features. 

As in DS we shall make one choice of Prandtl number (o=0.72), 

Mach number (M«=3.8) and ratio of specific heats (y=1.4). 


3.1. Cooled Wall Results 

We begin by considering the effect that wall cooling has on the 
inflexion points. In DS a condition was derived for the existence of 
inviscid, neutral, non-axi symme t ri c modes in the boundary layer on an 
axisymmetric body, which was termed the "triply generalised inflexion 
condition". However this condition does involve the azimuthal and 
streamwise wavenumbers a and n, and so it is difficult to forecast, 
prior to any numerical investigation, the existence of neutral stability 
points of this kind. However in the case of axisymmetric disturbances, 
this is no longer the case, since the condition reduces to 



w o r ( r ) 

rT 0 (r) 



0 , 


(3.1) 


as determined by Duck (12). 

Figure 2a shows the axial variation of radial position of the 
generalised inflexion points for the temperatures shown. As in the 
insulated cylinder case (12). the graphs display two prominent features: 
(i) the inflexion points occur in pairs and (ii) there exists a 
critical value of C. downstream of which no such points exist. The point 
; - 0 corresponds to the tip of the cylinder and as such corresponds to 
the planar case as studied by Mack (2). We find that as the surface 
of the cylinder is cooled, the lower inflexion point lifts up off the 
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cylinder surface. For sufficient cooling the lower point coalesces with 
the upper inflexion point and further cooling results in the complete 
disappearance of the inflexion points. Therefore for a given ^-station 
there exists a critical wall temperature below which no inflexion points 
exist. In light of our earlier comments, we do have an additional condition 
for the existence of neutral subsonic disturbances, namely that 

1 - L <c < ' + K, < 3 - 2 > 

This has direct implications on the first mode of instability as this 
generally requires the presence of an inflexion point in the profile 
sa t i s fy ing (3.2). 

Figure 2b shows the axial variation of wg(qj) for the displayed 

wall temperatures. From stability theory, unstable subsonic modes only 

exist if a generalised inflexion point satisfying (3.2) occurs somewhere 

in the boundary layer. Examination of the curve for T w = 3.0 reveals that 

subsonic generalised inflexion points only occur for 0 S ( < 0.0795; 

consequently the mode has completely disappeared before the generalised 

inflexion points have merged. For a wall temperature of T w = 2.0 the 

generalised inflexion points are all now supersonic in nature (c < 1 - 1 ), 

M» 

implying for this and all cooler wall temperatures the eradication of Mode I 
instabi 1 i t ies . 

We now turn our attention to the eigenvalue problem for both 
axisymmetric and non-axi symme t r i c disturbances. We shall present only 
results for unstable modes and all plots are for Cj (where c = c r + icj). 
The first set of results we present corresponds to the tip of the cylinder 
and as such are comparable to the planar results as obtained by Mack (2, 10, 
11). Figure 2c displays distributions for the first mode of 
instability. We observe that all the modes originate as sonic (i.e. with 
c = 1 * j|j-) neutral modes and terminate as subsonic, generalised 
inflexional modes. We note that as we cool the surface of the cylinder 


this mode undergoes stabilization, until, with sufficient cooling it becomes 
completely stable, thus verifying Mack's observations (and those of Lees 
(3), Van Driest (4) and Dunn and Lin (5)). From our inflexion point 
results, this is to be expected, as we observed that for cool enough 
wall conditions, w o(*li) drops below 1 - and we no longer have the 
necessary conditions for a subsonic generalised inflexional mode. 

Figure 2d shows the second mode of instability at the lip of the 
cylinder. All the modes originate as subsonic generalised inflexional 
modes with the special case c( = c r = w o(tl i ) > = 1 • which therefore 
corresponds to a generalised inflexion point in the frecstream, and 
terminate as neutral modes. Depending on whether or not a subsonic 
generalised inflexional mode for the given wall conditions exists, this 
neutral mode is either subsonic inflexional or supersonic, in nature, 
and the mode may continue as a decaying mode (Cj<0). Examination of our 
results reveals that as we cool the cylinder surface, the maximum value 
of c i increases slightly and then decreases but the important product 
acj , actually increases with wall cooling. Thus we deduce that wall 
cooling destabilises the second mode of instability, in line with Mack's 
(2, 10, 11) observations. 

The next set of results we present corresponds to a relatively small 
distance from the tip along the cylinder in the axial direction, at the 
location C = 0.05, and for an axisymmetric mode (ie, n=0) . Figure 2e 
displays the Mode I distributions and we note that with sufficient cooling, 
it is again possible to completely stabilize this mode. We observe that the 
neutral mode at which this instability originates now occurs for a 
slightly greater than zero and it is very slightly supersonic in nature. 

When compared with the corresponding planar results we note that curvature 
(as noted by Duck (12) and DS) has a profoundly stabilizing effect. Even 
though this station is only a relatively short distance along the cylinder 
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f rom the tip, curvature has reduced the value of maximum Cj for the 
T w = 3.0 curve by a factor of 4, while for the T w = 2.8 curve this 
factor is almost 5. Curvature results in the mode requiring less cooling 
to completely stabilize it. Figure 2f displays the axisymmetric mode II 
instability at this axial location. Again we note that curvature has had a 
stabilizing effect on the instability, but cooling causes the mode to become 
more unstable, in line with the planar results described previously. 

The next set of results we present corresponds to the £=0.05 location, 
for azimuthal wavenumber n=l . Figures 2g, 2h display the mode I and II 
instabilities, respectively. We observe that in this case the mode I 
instability is substantially more unstable than the axisymmetric case. We 
find again, however, with sufficient cooling we can completely stabilize 
this mode, although the mode does persist for cooler wall conditions. The 
mode II instability has the same qualitative features as the axisymmetric 
case, although it is slightly more stable than the axisymmetric case. 

In DS it was noted that near the cone tip, for non-axi symmet ric 
disturbance terms, a third mode of instability is seen to develop. This new 
mode originates as a neutral mode at a=0, with Cj*0, terminates as a 
supersonic neutral mode and was termed in DS, mode IA. As has already 
been noted, the work of DS was for a cone with adiabatic wall conditions 
and it was found that this new mode had already amalgamated with mode I 
instability at the £=0.05 location for n=l . We find however that for a 
cooled cylinder at this axial location and n=I , the mode IA instability 
is still distinct (Figure 2i) over the range of T w shown; a possible 
partial explanation is provided in Section 6. We observe that wall cooling 
causes mode IA instability to become less unstable and with sufficient 
wall cooling we can completely stabilize it. Comparing with the mode I 
instability we note that the maximum value of cj for the mode IA 
instability is an order of magnitude larger, but due to the smallness of a. 
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Ihe mode I growth rates are generally larger. One further observation, 
though, is that mode IA instability persists for cooler wall temperatures 
than does the mode I instability. 

We now consider the situation for an azimuthal wavenumber of n=3, at 
this axial location. We observe that for this case the mode I and IA 
instabilities have now amalgamated. The new combined mode originates as a 
neutral mode (but with Cj*0 for cc=0) and terminates as a subsonic 
generalised inflexional mode (Figure 2j). Again we note that sufficient 
wall cooling can completely stabilize this mode, although it is found that 
this mode persists for cooler wall temperatures as n increases. Comparing 
with the n=l results we find that the increase in n has also caused the 
mode to become slightly less unstable. The mode II instability (Figure 
2k), again, has the same qualitative features, but has also undergone 
stabilization with the increase in n. However, it is found that cooling 
has the more dominant destabilizing effect here. 

The next set of results we present correspond to an azimuthal wavenumber 
of n=5, at this axial location. We observe the same qualitative features 
as the n=3 results for both modes (Figures 21, 2m), although we note that 
we have had a more significant stabilizing effect due to the increase in n. 
The mode I instability is now completely stabilized for higher wall 
temperatures, while again, it is observed that cooling causes an even more 
marked destabilizing effect for the mode II instability. 

We find that as the azimuthal wavenumber n is further increased, both 
the mode I and II instabilities undergo additional stabilization, 
although cooling maintains a destabilizing effect on Mode II. For large 
n it is found that Mode I instability is completely stabilized while the 
mode II instability persists, but with much diminished growth rates. In 
the next section we shall consider the asymptotic structure of the 
disturbance equations in this limit. 
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The final set of results we shall present in this sub-scction corresponds 
to the axial location £=0.5 and for an azimuthal wavenumber of n=l . We 
find that at this distance along the cylinder the mode I and IA 
instabilities have now combined (Figure 2n). Comparison with the n=l 
results at the £=0.05 location, we note that growth rales have been 
reduced due to the stabilizing effect of curvature. It is found, however, 
that the combined mode prevails for cooler wail conditions, but again is 
completely stabilized with sufficient wall cooling. Figure 2o displays the 
mode II instability. Curvature has resulted in the growth rates being 
reduced, but we observe that cooling has a more profound destabilizing 
effect here. 

As we move further along the cylinder in the axial direction, the mode 
I and II instabilities undergo further growth-rate reductions. For a 
given wall temperature, T w , azimuthal wavenumber n, axial wavenumber a, 
there exists a critical value of £, beyond which no triply generalised 
inflexion points (see DS for details) occur. Consequently, we expect the 
mode I instability to have disappeared for axial distances larger than 
this critical £, which is bourne out by our numerical results. It is 
found however that due to the special nature of the inflexional point within 
the freestream, the mode II instability still originates with c=l , even 
for axial distances greater than this critical £, although the mode has 
much reduced growth rates. In the light of this, in Section 5 we shall 
consider the form of the disturbance equations in the limit of large £. 

3.2 Heated Vail Results 

We begin by considering the effect that wall heating has on the 
inflexion points.. We restrict our study to the case of axisymmetric 
disturbances, and consider the effect wall heating has on condition (3.1). 
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Figure 3a shows the axial variation of (radial) position of the 
generalised inflexion points for the temperatures shown. We observe again 
ihe features seen in Figure 2a. Close to the cylinder tip, however, we find 
that for a small axial distance measured from the lip, we no longer have any 
lower inflexional points. As we heat the surface of the cylinder further, 
this axial distance is found to increase. It is also observed that wall 
heating causes the critical value of beyond which no inflexional points 

exist, to increase. For T w =4.5, the critical value of C is about 
0.216, while for T w =6.0, we have had a substantial increase to a value of 
£ a 0.423. This will have direct implications on the first mode of 
instability which we expect to persist for larger axial distances as the 
wall is heated. These effects are in many ways to be expected, being the 
converse of the cooling observations described earlier. 

Figure 3b shows the axial variation of WQ(Tli) for the displayed wall 
temperatures. The most marked feature of these curves is that as we heat 
the cylinder surface, the lower inflexional point becomes subsonic past a 
critical value of C, which is temperature dependent, i.e. for axial 
distances greater than this critical value of C but upstream of the station 
beyond which no inflexional points occur, both generalised inflexional 
points are now subsonic in nature. The first wall temperature for which we 
observe a lower, subsonic inflexional point is for a wall temperature of 
about T w =4.5. It is found, however, that the critical value of C here, is 
very close to the stations where the generalised inflexional points coalesce 
For T w =5.0 we observe that for the range 0.2635 < C < 0.2720 we have two 
subsonic inflexional modes, while for the hotter wall temperature of T w =6.0 
we have the larger range 0.363 < £ < 0.423 for which both generalised 
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inf Icxional points arc subsonic. In these £-ranges we expect the appearance 
of two subsonic generalised inflexional modes to have a significant effect 
on the physics of the problem. 

We now present eigenvalue results for axisymmetric disturbances only. 

As in the case of cooled wall conditions we focus attention on unstable 
modes. We begin by considering the effect wall heating has on the 
mode I and II instabilities for a ^-station close to the cylinder tip 
(£=0.05) and consequently the lower inflexional point is still supersonic 
in nature. Figure 3c shows the mode I instability for the temperatures 
shown. We observe that all the modes originate as neutral modes at a value 
of a slightly greater than zero, and are very slightly supersonic in 
nature. We find as we heat the wall this neutral mode tends towards the 
sonic value. All the modes terminate as subsonic generalised inflexional 
modes, continuing on as stable modes ( c j <0 ) for larger values of a. 

These observations are similar to the results obtained for the axisymmetric 
cooled wall case for this ^-station. We find, as we would expect, heating 
the surface of the cylinder causes the mode I instability to become more 
unstable - converse to the effect of cooling on this mode. 

Figure 3d displays the mode II instabilities for the £=0.05 station, 
at the temperatures shown. We find all the modes originate as subsonic 
generalised inflexional modes, rise to a maximum and terminate as subsonic 
generalised inflexional modes (which then continue in all the cases 
presented as stable modes). Heating the cylinder wall causes the mode II 
instability to become less unstable and the numerical evidence suggests that 
with sufficient heating we can completely stabilize this mode. 

We new consider the effect that the lower inflexional point becoming 
subsonic has on the mode I and II instabilities. Figure 3e displays 
the mode I instability for a wall temperature of T w =5.0 and for the 
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£-stations as marked. For the £=0.26 station (where the lower inflexional 
point is still supersonic) the mode originates as a very slightly supersonic 
mode and terminates as a subsonic generalised inflexional mode. For the 
£=0.264 station the lower inflexional point has now become subsonic in 
nature, and we find the mode I instability now originates as this lower 
subsonic generalised inflexional mode. Consequently the value of a for 
the neutral mode has increased correspondingly. As before, the mode I 
instability terminates as the upper generalised inflexional mode which is of 
course subsonic, as well. From the inflexional point curves we know that as 
we move upstream the inflexional points move closer together, eventually 
coalescing and this is reflected in the new form of the mode I 
instabilities. For the £=0.27 stat ion the mode I instabi 1 i ty occurs 
over a much smaller a-range and the growth rates are greatly reduced. 

Figure 3f displays the mode I instability for a wall temperature of 
T w =6.0 and the displayed £-stations. Again we note that as the lower 
inflexional point becomes subsonic the neutral point at which the 
instability originates transforms from being very slightly supersonic in 
nature, to this inflexional mode. Movement upstream causes the a-ranges 
and growth rates to be diminished, but the reduction is less marked (in 
comparison with the T w =5.0 results) due to the destabilizing effect 
brought on by wall heating. 

The appearance of a second subsonic generalised inflexional mode is 
found to have no effect on the mode II instability as it always terminates 
as the upper generalised inflexional mode. It is found, however, that a 
third mode of instability exists, originating as the lower generalised 
inflexional mode and terminating as a slightly supersonic neutral mode. 

This new mode, which we shall term Mode IIA, occurs for values of a 
greater than the value of a for which the mode II instability terminates. 

It appears that the mode II instability continues as a stable mode and 
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thcn becomes unstable again at the lower generalised inflexional point. 

The growth rates of the mode 1 1 A are found to be very small. For a wall 
temperature of T w =5.0 and the station £=0.27 the growth rates are in the 
order 10“ll t while for T w =6.0 and £=0.364 the growth rates are of the 
order of 10" 13 - 10" 14 . 

We now move on to consider the form of the disturbance equations in 
the limit of large azimuthal wavenumbers. 


§4. Disturbance Equations for Large n 

In this section we consider the form of the disturbance equations in the 
asymptotic limit of large n, guided by our numerical results. 


4.1 Formulation of the Problem 

Consider the pressure disturbance equation as derived in DS 


4- [ ] + r _ W()Tl] J&ln 

a 2 dn L w 0 -c J L 1+TiC ( * lJ a 2 (w 0 -c) 


= 4»P, 


( 4 . 1 ) 


whe re 


* = T ° C 1 + 1 " M - {w 0- c > 2 - (4-2) 

In the limit of large n our numerics suggest that the corresponding 
streamwise wavenumbers for the instability also increase, and that 
4 >(ti= 0) — * 0. Equation (4.2) suggests we must have 


a = an, a = 0(1). (4.3) 

A (sensible) balancing of terms in equation (4.1) suggests asymptotic 
expansions of the form 

c = cq + n* 2/3 cj + 0(n' 4/3 ) + ... , 

<J> = Op + n* 2 ^ 3 4>j + 0(n* 4 ^ 3 ) + ... . 


(4.4) 
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The leading order pressure equation has the form 

,2 


*Vi “ T5 p " °' 


(4.5) 


where 


«o' T o [ 1 * a z |i ! nCI 2 ] - M - <" 0 -‘ 0 > 2 - 

This has a solution of the (bounded) WKB type (assuming 
> 0 

, • . * 

~W t To 1 

We can re-write the expansion for <t> in the neighbourhood of q 


(4.6) 


p - -*Sr { ' TST f [ - *0 ]' "1 )■ 


(4.7) 


= 0 as 


<J> = d> 0 (0) + n' 2/3 [ <t> 0 (0) q + 4>, (0) ] + 0(n' 4/3 ), 


(4.8) 


where 

4*o(0) = T w [ 1 + -Aj ] — M_ c n , 

*1(0) = - 2 c 0 c, M 2 , 
and for insulated wall conditions 


(4.9) 

(4.10) 


= _ 2 A 3 la + 2 M 2 w' 


<yo) = - z -^r r JL + 2 M - w o (0) c o- 

while for heated or cooled wall conditions we have 


(4.11) 


V°> - w°> c ' * % 1 - 4^ + 2 M - V°> V 

where T w represents the wall temperature. 

If 0g(0) = 0, then we must have 


(4.12) 


c 0 ■ 


■ t i [ i + i ]* 

w l q2 j 


TC 


(4.13) 


(which is clearly real). We now seek to determine the first order 
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correction term to c, namely cj, For this we must study the wall layer, 
which by a balancing of terms of (4.2) suggests is given by 

- -2/3 

T 1 = ti n , tl = 0(1) . (4.14) 

The (bounded) solution for P in the wall layer is given by 

P = P 0 Ai(ii), (4.15) 

where Pq is a constant and q is given by 


Tl = 



-i 1 / 3 _ ■ _ 

| [ * 0 (0)n + <t>,( 0 ) ] . 


(4.16) 


The boundary condition at the wall is 


tl=0 " 0t 


ri i In 


which we transform to rf-space, giving the equation 
P- [X l/3 ij.B] =0. 


(4.17) 


(4.18) 


where 


x. — t- 


T w p 0 (0)] 


7 . B = OjCO) 


(4.19) 


Now since the solution in this region to the pressure is given by Airy's 
function (4.15) we have that 


Ai q [n = BX ] = 0. 

Transforming back to q* s P ac e we have the result 

1 1/3 


(4.20) 


Yn - 


a 2 


T w [%( 0 )] 


[ - *1(0) ], 


(4.21) 


where y n (where n » 1,2,...) represent the solutions of the equation 


Ai'(- y n ) = 0. 


(4.22) 
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Substituting either (4.11) or (4.12) anti equation (4.10) into equation 
(4.21) yields the first order correction term lor C\ namely 


C 1 


t 1/3 
1 w 


( 2 p: 




c„ - 


l\v V 


2 a 


TTT 




n 


!/3 


Yn- 


(4.23) 


for insulated walls, whilst lor lie a led/eoolcd wall conditions we have 

2 pi v°» c n - ^ 3 * v 1 ” 

2 fi2/3 c 

oo l) 

where cj is real for both cases. 

Since y n is a solution ol equation (4.22). where y n > 0, then we 
have an infinite number of discrete, real possible values for y n as the 
derivative of the Airy function has an ini ini to number of discrete roots, 
confined to the negative real axis. This suggests that there arc an 
infinite number of modes. 

We shall now compare these results with numerically determined results 
for large values of n. 



,2/3 


Yn- (4.24) 


4.2 Numerical Results 

All the results presented in this section arc for a frccstrcam Mach 
number of 3.8 and at the point £ = 0.2 along the cylinder. We only 
present results for adiabatic wall conditions. 

Firstly let us consider the asymptotic expansion for c in the limit 
of large n. The leading order term in the c expansion cq. is given by 
equation (4.13) and using the numerically determined values 


T w = 3.379, 
a = 0.1525. 


(4.25) 


we have 


c 0 = 0.7978. 


(4.26) 
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The first order correction term, which for adiabatic wall conditions 
is given by equation (4.23), is computed to have the value 

c | = 0.3698 y n , (4.27) 

where we have used wq^O) = 0.1904 and y n are the zeros of the 
derivative of the Airy function. The first six possible values of y n are 
determined from tables (Abramowitz & Stegun (16)) and the corresponding 
values of cj are shown in table I 


C 1 


0.3698 

(I) 

1.2092 

(II) 

1 . 7750 

(HI) 

2.5146 

(IV) 

2.9214 

(V) 

3.3282 

(VI) 


TABLE I: VALUES OF cj 

-2/3 

Figure 4a shows a plot of c( = cq + n cj), as determined 
asymptotically, against n for the different values of the first order 
correction term cj, where the numbering refers to the numbering of the 
correction terms in Table I. (It should be noted that only n integer has 
physical significance, although Figs. 4a and 4b show c as a continuous 
function of n.) From hereon we shall refer to the different values of 
c as order I to VI inclusively corresponding to the numbering convention of 
the correction terms in Table I. 

Now as observed above, our asymptotic analysis suggests the existence 
of an infinite, discrete number of possible values for c. When we searched 
for the eigenvalues numerically, for large n, we determined that there 
were indeed many modes. Figure 4b displays two plots of c r against n 
for order I and order V correction terms. Graph (1) in each case represents 
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the asymptotic curve and graph (2) is the numerically determined curve. It 
should be noted that in this range of n and a, iCji « 1 
(cj ~ 10- 10 — 10 ' 1 2 ) , comparable to the machine accuracy of our 
computations. From our two sets of plots we note that we have good 
agreement between the numerical solutions and asymptotic theory for large 
n . 

Turning our attention now to the form of the pressure disturbance terms, 
as obtained numerically, we find that they do indeed follow the pattern 
predicated by our asymptotic theory, being initially oscillatory in the 
Airy solution region and then decaying to zero in the far field. It is 
also observed that increasing the order of the correction term has the 
effect of increasing the number of zeros of the eigensolut ion. Figure 4c 
displays the distributions of Real { P } for n = 40 corresponding to 
the orders as shown. 

Examining the c expansion again, we have determined that both cq 
and Cj are real and therefore the leading order imaginary term. Cj, 
is at most 0(n‘ 4/3 ). This means that the leading order term in the growth 
rate (a cj) is 0(n' 1/3 ) at most. Therefore actual growth rates will 
decrease as n — ♦ ~ which is confirmed to be true by our numerical 
observat ions . 

We now turn our attention to the form of the disturbance equations in 
the limit C — * 00 ■ 


§5 Disturbance Equations for Large £ 

In this section we consider the form of the disturbance equations in 
the far downstream region, guided by our numerical observations. 

5,1 Formulation of the Problem 

Consider the pressure disturbance equation, as presented in the previous 
section, but written in terms of 'r' rather than f q', 
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^Q -Zi. . C . 2 = d,p ( 


( w o* c ) 


(5.1) 


whe re 


* = T 0 [ 1 + ^77 ] - M f ( w 0' c > 2 • 


(5.2) 


and 


r = I + 1) C- (5.3) 

In the limit ^ we assume a scale of a of the form 

a=^C. (5.4) 

where a* is to be determined. 

Guided by Duck s (12) work for the form of the basic flow in the 
farfield we define the (small) parameter 

« ■ G log Z) * 1 = (log O' 1 - (5.5) 

In the slow moving viscous region close to the wall (namely the r=0(l) 
lengthscale) we expect asymptotic expansions of the form 

c = c 0 + e cj + 0(€ 2 ), 

A A 

4> = <J>o + e 4>i + 0(e 2 ) , 

T 0 = T w + « Tj + 0(e 2 ) , 

w 0 = < w 0 + 0(e 2 ). (5.6) 

Now, our numerical observations strongly suggest that <J>g — * o as £ 
becomes large in the r = 0(1) region (which may be confirmed a posteriori) 
which means that = 0(e). Examination of the <& expression reveals that 
this is only possible in general if 

= 0(e-i). (5.7) 

Therefore we can redefine the scale (5.4) to have the form 

« = « C (log O*. (5.8) 
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whcre 

a = 0(1). 

To leading order, equation (5.2) reduces to 

*o = T w - "1 V (5 9) 

but in view of our comments above that <J>q — * 0, as C, —* for r = 0(1). 
then we have 

Jo-fe 1 - (5 "» 

which means cq is real. 

At the first order in « equation (5.2) has the form 

it = T\ + nil* - M* (2 c 0 cj - 2 c 0 w 0 ), (5.11) 

while the 0(«) correction to the pressure equation (5.2) is given by 

Vrr J+P, * [2 l 0 »l (i, - I 0 ) -T, - gS] P-0. (5.12) 

We now transform equation (5.12) using the same transformation as 
applied by Duck (12) for the basic flow. Firstly we employ the transform 

7 = lnr, (5.13) 

giving 

Ife- 2 'P„* [ 2 ; 0 (2l ■ ^o) -T, -2^e- 2r ] P-0. <5.14) 

Secondly we use the transform 

(515) 

J 0 w 

but since T w is constant with respect to r we can simplify (5.15) to 
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f = IT T w , (5.16) 

where we take the constant of integration to be zero. Equation (5.14) can 
be re -wr i t ten 

P RR + { ® 2 T w [ 2 c 0 (ci - K) - T,] e 2 * T » - n 2 T 2 } P = 0, (5.17) 

where we have made use of the result obtained by Duck (12) for the basic 
flow 


w 0 = (5.18) 

Equation (5.17) is solved numerically to obtain a value for cj subject to 
the condition at the wall 



(5.19) 


and that P is bounded in the far-field. The second condition is obtained 
by taking ~ limit of (5.17), i.e. 


2K T„ 


where 


Pj^ - p 2 K e'- 1 ' P = 0, 
ft 2 = 2 c 0 M 2 a 2 T w . 


(5.20) 

(5.21) 


We find, to leading order, this equation has a decaying solution of the 

form 


P - |T* K ' 1/4 exp [ - if- e R T *R i _ \ IT T w ]. (5.2 

5.2 Numerical Results 

All the results presented in this section correspond to a free-stream 
Mach number of 3.8, and azimuthal wavenumber n=l . 
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From our numerical observations we obtain a value for the leading order 
term in the c expansion (5.6) of the form 

cq = 0.4617922. (5.23) 

Using a fourth order Runge-Kutta scheme equation, (5.17) was solved 
subject to conditions (5.19) and (5.22) to determine the eigenvalues cj. We 
find that for a given value of ft there appears to be a large number of 
discrete, real values for cj. Figure 5a displays a plot of cj against 
ft corresponding to the first five modes, as shown. 

As observed above, our asymptotic analysis implies the existence of 
a large number of discrete possible values for c \ , which in turn implies 
the existence of a large number of discrete values for c r . When we 
searched for the eigenvalues by solving the full system of equations 
numerically, for large £, we determined that there were indeed many modes 
and we managed to identify the first five modes. Figure 5b displays a 
comparison between the asymptotically determined value of c r and the 
numerically determined value of c r against C corresponding to the first 
mode. We have relatively good agreement, since the error term in the 
asymptotic theory is 0(«2), which is quite large. Therefore the numerical 
results seem to confirm our asymptotic theory. 

The asymptotic theory presented above tells us nothing about Cj and 
consequently reveals no information about the growth rate etc,; such an 
investigation would require a prohibitive amount of algebra. 

We now move on to consider the form of the disturbance equations in the 
limit of small (large) £ for cooled/heated wall conditions on a cone. 
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§6 Disturbance Equations for small £ 

In this section we consider the effect wall heating/cooling has on the 
asymptotic theory developed by DS for the limit £ —t » 0 (and also with 
some simple modification the limit C — » “ in the case of a cone) for the 
description of mode IA. It should be noted that the analysis in this 
section is also applicable for a cone, the circular cylinder case being 
retrieved when X = 0. We find that when we impose heated/cooled wall 
conditions, instead of adiabatic conditions as treated by DS, the asymptotic 
theory is rather different. Here we consider the case of a — * 0, £ — » 0, 
In DS it was shown that the disturbance equation for <p in the limit 
ct — * 0 (assuming n * 0) reduces to 



In the limit £ — * 0, we assume expansions for the q = 0(1) scale of 
the form 


the form 

C = Ccj + C 2 C2 + C 3 C3 + . . . , (6.2) 
9 = 90^) + C9l(Tl) + C 2 92(H) + C 3 93(H) + ... , (6.3) 
w 0 = W 00( T 1) + *>0lOl) + C 2 w 0 2(n) + C 3 W03Ol) + ••• . (6.4) 
T 0 = TqoOI) + C T 0 iOl) + C 2 Tq2( t 1) + C 3 To3(H) + ... . (6.5) 


where wqo(ti) and TgoC 1 !) represent the planar values of the velocity and 
temperature profiles respectively, and w 0 i(q) and Tqi ( tl) etc. correspond 
to the perturbations to the basic flow caused through curvature. A more 
complete expansion for <p is given in DS. 

To leading order, we have for q = 0(1) that 


90 = A 0 w 00( | l>- 


( 6 . 6 ) 
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whcre Ag represents an arbitrary amplitude parameter. 

At the next order in £ we obtain the equation for <pi 

w 00 9 1 q “ A 0 w 00q *1 + A 0 w 00 2 “ w 00q 91 
+ Aq w 01 w QOn - Aq w 0]ll w 00 = kj Tqo • 

where we have used equation (6.6) and k| is a constant. Setting q = 0 
in equation (6.7) and assuming i J ^=0 = 0 wc havc 

- A 0 w 00q ( T l = °) = k l T 00 (n=0) . (6.8) 

The boundary conditions as q — » <*>. must be compatible with the form 
of system (6.1) in the far-field (as presented in DS) , together with (6.6). 
Defining 

7 = 1 + q C = OO). ( 6 - 9 ) 


then we must have an outer solution of the form 


9 


out 

0 


= A f 


n- 1 


( 6 . 10 ) 


where 


A = A 0 + C A 1 + ■ • • • 
and so we must also have 

^lqlq , ~ ~ 9of | f = i = *( n+ l) A 0- 


( 6 . 11 ) 


( 6 . 12 ) 


Substituting this into (6.7) yields 

(6.13) 

In order to obtain complex values of c in our asymptotic analysis we 


- n Tqo (n=Q) 
1 w 00q (n=o) 


must consider higher orders in £. 

At the next order of £ we have the system 
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w 00 92q + 2 * w 00 <POn + 2t 1 w OO 9 1 ti + w 02 9oq 

- w 02t| 90 + W 01 <Plri + I 2 w OO ^Otj “ cj <P|r| 

+ 2 W 01 ^ 90q + W 01 90 “ w 01ri 9l - 2, 1 w 0 1 rj 90 

- 2r\ ci <POr) - 90r| + w 00 9l + *1 w 00 90 

- cj <P0 “ w 00q 92 “ 2 ^ w 00q 90 ~ 2r \ 9l w 00q 

T) 

- t 1 2 w 00n 90 = k 2 T 00 + ° 2 T 00 f ^° t ^ Q dT l 

J Q ‘00 

+ *1 Tqj . (6.14) 

We now determine the leading order imaginary component of the complex 
wavespeed c. Since (6.14) contains only real coefficients, any 
imaginaries must, of necessity, only arise at a critical point, where 
c = wqq. Since c = 0(0* this must occur when i\ = 0(0. We therefore 
consider a thin layer relative to the T| = 0(1) scale, namely 

y = ti/C = 0(i). (6.15) 

On this scale, the expansion for <p is expected to develop as 

9 = C ^o (7) + C 2 ^i (y) + C 3 $2 (y) + • • • • (6.16) 

where the <l>j are expected to be normalised in such a way as to be 
generally 0(1) quantities. 

The leading order equation has the form 

*0 ' A 0 y w 0(hl (q = 0) (6.17) 

(where Aq was introduced in equation (6.6)). 

At the next order in £ we obtain the system 

[ w 00q(°) 7 * cjJ 3>iy - w 00n (0) = kj 

+ i V 2 *00qq(0) A 0 w 00n (0) - T 00n (0) y. (6.18) 

where i(] is a constant term. This system cannot be simplified any 
further, as in DS, since the condition Woo^(q=0) = 0, no longer holds 
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in the present (i.e. hcatcd/coolcd wall) ease. 

As in DS , we assume Ag and kj to be real quantities, as we arc 
at liberty to do so, and now consider (where i denotes an imaginary 

component). This quantity is triggered by the +itt jump in the logarithm 
across the critical layer. We can re-write (6.18) in the form 

[ w OOt|(°) 7 - ci ] *iy - $1 w 00r|(0) = R - (6 ' 19 


where 

R = ki + i y 2 w 00 t|T|( 0 ) a 0 w 00 t|(°) “ A ° C { 0 o(8)“ )t ° 0t 1 (0) 

( 6 . 20 ) 


and find 


$1 = ( w 00q(°) y ‘ 


f 


R dy 

o [ w ooti(°) y- c i^ 


( 6 . 21 ) 


Evaluating this integral, taking just the imaginaries together with the 
1 i mi t as y — » « yields 


~ K A 0 C! 


( w oon(°) y-cil J 


w 00nn(01 T QQn (0 J 

w 00 2 ti(°) t 00(°) w 00ti(°) 


) 


(6.22) 


= B' [W 00t 1 (0) y - C!l. 

where B< = s A 0 c, { } . (6.23) 

l W00 2 q(°) T 00(0) w 00q( 0 ) J 

Equation (6.22) now provides a lower boundary condition for systems (6.7) 
and (6.14). 

Assuming (6.7) contains real coefficients only (c j is assumed real, 
an assumption that may be justified a posteriori) we have 


<P 


i 

1 



W 00(B). 


(6.24) 


where B* is defined in (6.23). 
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Considcring the imaginary part of system (6.14) in the limit q — * « 
we obtain 



+ v 

1 



(6.25) 


By the form of <p in the far- field (as presented in DS) we require 


<P 


i 

2n 


- - (n+l ) <p! I 

« 1 I f| ) 00 

- - (n+l) B' , 


and therefore 

= - B'n. 


(6.26) 


(6.27) 


We now consider the imaginary part of (6.14) evaluated at q = 0 


• ‘l finite - J 2 % U - w ooi<°> tlU - k 2 T oo<°>- 

Using the results 


we obtain 




w 00q(°) • 



= - B' c, , 


c\ = n 2 g. T QQ^ =0 > / w QQrm( T l = Q) _ iQQnlllfg) l. 
W 00 3 t 1 <tl=0) ' w 00tl(^l=0) Too(U=0) j 


(6.29) 


(6.30) 


Using the governing equations of motion we can simplify (6.30) to the form 

cl = n 2 — T QQ 2 j 0) *QQrni(°), (6 3] 

wg On (0) 

Comparing this result with the adiabatic result, as obtained in DS, 
we note that the first imaginary term in c is an order in £ larger for 
heated/cooled wall conditions, implying larger growth rates in the present 
case. The ratio of the leading-order imaginary terms has the form 
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Heated/Cooled Cj _ 2 w OOnn(Q) . (6.32) 

Adiabatic Cj C w 00r)(0) ( 1+ 2o (Y-l) n l 

We find, however, that altering the wall conditions from adiabatic 
conditions as treated by DS. to hcated/cooled wall conditions, has little 
significant effect on the asymptotic analysis for the mode I lower neutral 
point for which a = 0(^1). 


§7. Conclusions 

In this paper we have investigated the supersonic boundary-layer flow, 
and the inviscid stability thereof, over axi symmet ric bodies (with 
particular emphasis on a long thin, straight circular cylinder) subject to 
heated or cooled wall conditions. 

The effect of wall cooling is, generally, seen to reduce the importance 
of the ''first mode'' of inviscid instability, while the amplification rates 
of the ’’second mode'' of inviscid instability are increased. Therefore we 
have agreement with the effect wall cooling has on planar boundary layers 
(Mack (11), for example). 

The converse effect is observed with wall heating. The amplification 
rates of the "first mode" of instability increases, while wall-heating 
causes the "second mode" of instability to be stabilized. Wall-heating 
may also cause the formation of a second subsonic generalised inflexional 
mode, which results in the appearance of an additional mode of instability, 
not found in adiabatic or cooled-wall studies. 

Asymptotic studies for large azimuthal wavenumbers and far downstream 
distances reveal that the real part of the eigenvalue c is non-unique 
in both limits, suggesting that there exists an infinite number of discrete 
possible values for the real part of c in both these limits, although the 
corresponding values of the imaginary part of c are exceedingly small. 
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